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ABSTRACT 

We present 350 /im observations of 15 Chapman et al. submillimeter galaxies (SMGs) with radio counterparts 
and optical redshifts. We detect 12 and obtain sensitive upper Hmits for three, providing direct, precise measure- 
ments of their far-infrared luminosities and characteristic dust temperatures. With these, we verify the linear 
radio-far-infrared correlation at redshifts of z ~ 1-3 and luminosities of lO^^-lO^^L©, with a power-law index 
of 1.02 ±0.12 and rms scatter of 0.12dex. However, either the correlation constant q or the dust emissivity 
index (3 is lower than measured locally. The best fitting ^ :^ 2. 14 is consistent with SMGs being predominantly 
starbust galaxies, without significant AGN contribution, at far-infrared wavelengths. Gas-to-dust mass ratios 
are estimated at 54;^}^ (/^850/im/0.15 m^ kg~^) , depending on the absoption efficiency Kj,, with intrinsic disper- 
sion ^ 40% around the mean value. Dust temperatures consistent with 34.6 ± 3K(/3/1.5)~^-^\ at z ~ 1.5-3.5, 
suggest that far-infrared photometric redshifts may be viable, and perhaps accurate to 10% < dz/il +z), for up 
to 80% of the SMG population in this range, if the above temperature characterizes the full range of SMGs. 
However, observed temperature evolution of T^ oc (1 +z) is also plausible and could result from selection ef- 
fects. From the observed luminosity-temperature {L-T) relation, L ex T^^^^^'^*^ , we derive scaling relations for 
dust mass versus dust temperature, and we identify expressions to inter-relate the observed quantities. These 
suggest that measurements at a single wavelength, in the far-infrared, submillimeter or radio wave bands, might 
constrain dust temperatures and far-infrared luminosities for most SMGs with redshifts at z ~ 0.5-4. 
Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: ISM — galaxies: photometry — 
galaxies: starburst — infrared: galaxies — submillimeter 



1. INTRODUCTION 

Ever since the first submillimeter- selected galaxy (SMG) 
samples debuted from SCUBA 850 /im surveys (Small et al. 
T997 [ lBarger"eran 11998: Hughes et al. 1998; Eales et alj 
1999h . the nature of the SMG population has been the fo- 
cus of attention for the galaxy formation community, be- 
cause the 850 /im selection is expected to pick similar sources 
almost independently of redshift (z ~ 1-8), due to a neg- 
ative ^-correction that essentially compensates for the loss 
of flux from increasing distance. This allows unbiased, 
luminosity-selected studies of galaxy formation. The hunger 
for information on these sources spurred a flurry of follow- 
up studies at all wavelengths, long and short. Since then 
many of the se sources have been identi fied at optical, UV 
(iBorvsetal .1120031: IChanman et aTl'2003^ "Webb et al.« l2003<) 
and radio wavelengths I Small et al. 2000; Ivison et al. 2002), 
providing accurate positions, which allowed optical redshift 
measurements (Chapman et al. 2003, 2005). As a result we 
now know that these massive galaxies, with redshifts dis- 
tributed around z — 2.3, are enshrouded with such quantities 
of dust that they often lie hidden at optical wavelengths, and 
therefore constitute a distinct population from the galaxies se- 
lected by optical surveys. 

More recently, longer w avelength subm illim eter surveys, 
at 1100 and at 1200 /im ("La urent et all 120051: lGreveetal.1 
12004), added to the pool of available information. However, 
the close proximity of the SCUBA, Bolocam, and MAMBO 
wavelengths on the Rayleigh- Jeans side of the spectral en- 
ergy distribution (SED) does not allow for an effective con- 
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straint on the thermal far-infrared spectral SEDs at the rele- 
vant redshifts. Nor do the latest results from the Spitzer Space 
Telescope provide powerful constraints, since at the shorter 
mid-infrared wavelengths the emission is dominated by poly- 
cyclic aromatic hydrocarbons (R\Hs) and a minority popu- 
lation of hot dust. For these reasons, the best estimates of 
the characteristic temperatures and the integrated luminosi- 
ties, to date, have relied on the assumption that t he local 
radio to far-infrared correlation (Helo u et al]ll985t ICondonI 
1992; Yunetal. 200 1) can be extended into the distant uni- 
verse. There_^ejiints that this may be appropriate ( Garrett] 
'2002UAppleton et ar.L2004) . but the assumption has remained 
largely unchecked. 

Shorter wavelength submillimeter measurements, sampling 
near the rest-frame peak of the emission, are thus essen- 
tial to provide firm constraints to the far-infrared SED. Here 
we present results at 350 /im, using the second-generation 
S ubmillimeter High Angular Resolution Camera (SHARC- 
2: lDowell et al.ll2 003) at the Caltech Submillimeter Observa- 
tory (CSO). From these we derive the first direct measures of 
dust temperatures and far-infrared luminosities for a sample 
of SMGs, testing the radio to far-infared correlation. We also 
attempt to constrain dust emission properties and investigate 
the implications of our findings for the viability of photomet- 
ric redshifts based on far-infrared and radio measurements. 
Finally, we present a range of useful scaling relations that may 
apply to the SMG population. 

2. SHARC-2 350 /im OBSERVATIONS 

We conducted follow-up observations of SCUBA 850 /xm 
detected sources with radio identifications and optical red- 
shifts (Chapman etal. 2003, 2005). Seven of the 15 tar- 
gets were hand-picked on the basis of their pre dicted bright 
350 /im fluxes (5'350/im > 45mJy from iChapman et al...2005.) . 
while the remaining were selected at random to partially com- 
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Fig. 1. — 30"x30" SHARC-2 350/im thumbnail images of the SMGs observed (sources 1-15 in Tables fTlsl shown in reading order (left to right, top to 
bottom). Images are centered on the radio positions, and displayed on the same flux scale for comparison. The 1274 FWHM effective beams of the optimally 
filtered (convolved) images are also indicated. Sources 3, 9, and 10 are not detected. 



pensate for any selection bias. 

The observations were carried out during eight separate ob- 
serving runs between 2002 November and 2005 April, in ex- 
cellent weather ijiis guz < 0.06), reaching la depths of 5- 
9 mJy in 2-4 hours of integration in 14 small fields (:^ 2.5 x 1 
arcmin^) around the targeted sources. Our scanning strat- 
egy was to modulate the telescope pointing with a small- 
amplitude (15"-20") nonconnecting Lissajous pattern within 
the limits of reasonable telescope acceleration (with typical 
periods of 10-20 s). This pattern was chosen to provide 
fast, two-dimensional, recurring but non-closed patterns with 
crossing paths — all of which are essential to allow the separa- 
tion of the superposed source, atmospheric, and instrumental 
signals. For the observations since 2003 February, we have 
taken advantage of the CSO's Dish Surface Optimization Sys- 
tem (Leong 2005) to improve beam shapes and efficiencies, 
and hence sensitivities, at all elevations. 

Pointing corrections were retroactively applied to the ob- 
servations by a careful fit to the pointing data taken during 
each run. In addition, our preliminary source identifications 
revealed a small, albeit significant (c^ ?>") systematic pointing 
shift in the negative R.A. direction. All maps were realigned 
accordingly. The reconstructed pointing is accurate to 3 '.'5 
rms, in good agreement with CSO specifications. 



The data were processed using the CR USH^ software pack- 
age, developed at Caltech (^Kovacs~2006), which models total 
power detector signals using an iterated sequence of maxi- 
mum likelihood estimators. The FWHM point- spread func- 
tions of the final maps (Fig.[l} are approximately 12'' after 
optimal filtering to search for beam-sized (c^ 9") features in 
order to yield maximal signal-to-noise ratios on point sources. 
This degradation of the image resolution by approximately 
a/2, from the nominal instrumental beam width of 8'f5, is be- 
cause the fitting for beam- sized features is equivalent to con- 
volving the image with the slightly wider effective beam of 
around 9" , which accounts for the smearing of the nominal 
beam by pointing and focus variations on the long integra- 
tions. 

Calibration was performed primarily against planets and as- 
teroids, when available, or using stable galactic continuum 
sources and Arp 220.^ Trends in aperture efficiency, espe- 
cially with elevation, were estimated and were taken into ac- 
count when calibrating our science targets. The systematic 
flux filtering effect of the aggressive reduction parameters 
used in CRUSH to subtract noise signals and reach maximal 
image depths, was carefully estimated for the case of point 

^ See http://www.submm.caltech.edu/ share/crush 

^ See http://www.submm.caltech.edu/ sharc/analysis/calibration.htm 
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TABLE 1 
Summary of Observations 



ID 



Name 



Offset^ 
(J2000.0) 
(arcsec) 



5(350 Aim) 
(mJy) 



5(850 /im) 
(mJy) 



5(1 100 /im) 
(mJy) 



5(1200 /im) 
(mJy) 



5(1. 4 GHz) 

(/iJy) 



1 


SMM J030227.73+000653.5 


+1.0, 


+0.8 


1.408 


42.2 ± 


9.8 


4.4 


± 


1.3 














217 


± 


9 


2 


SMM J105207.49+571904.0 


-1.0, 


+3.7 


2.689 


38.0 ± 


7.2 


6.2 


± 


1.6 








(0.4 


± 


0.8) 


277.8 


± 


11.9 


3b 


SMMJ105227.77+572218.2 






1.956 


(11.3 ± 


6.7) 


7.0 


± 


2.1 


5.1 


± 


1.3 b 


3.1 


± 


0.7 


40.4 


± 


9.4'^ 


4b 


SMM J105230.73+572209.5 


+3.3, 


-1.8 


2.611 


41.0 ± 


6.8 


11.0 


± 


2.6 


5.1 


± 


1.3 b 


2.9 


± 


0.7 


86.3 


± 


15.4 '^ 


5 


SMM J105238.30+572435.8 


+1.4, 


+2.5 


3.036 


40.5 ± 


6.5 


10.9 


± 


2.4 


4.8 


± 


1.3 


4.8 


± 


0.6 


61.0 ± 


22.0 '^ 


6 


SMM J123600. 15+621047.2 


-1.4, 


+2.0 


1.994 


22.3 ± 


6.3 


7.9 


± 


2.4 














131 


± 


10.6 


7 


SMM J123606.85+621021.4 


+6.7, 


+3.5 


2.509 


35.1 ± 


6.9 


11.6 


± 


3.5 














74.4 ± 


4.1 


8 


SMM J131201. 17+424208.1 






3.405 


21.1 ± 


7.7 


6.2 


± 


1.2 














49.1 


± 


6.0 


9 


SMMJ131212.69+424422.5 


+2.2, 


-4.2 


2.805 


(3.7 ± 


4.4) 


5.6 


± 


1.9 














102.6 


± 


7.4 


10 


SMM J131225.73+423941.4 






1.554 


(14.7 ± 


7.4) 


4.1 


± 


1.3 














752.5 


± 


4.2 


11 


SMM J163631.47+405546.9 


-1.3, 


-4.0 


2.283 


38.3 ± 


5.5 


6.3 


± 


1.9 








(1.1 


± 


0.7) 


99 


± 


23 


12 


SMM J163650.43+405737.5 


-1.0, 


-1.7 


2.378 


33.0 ± 


5.6 


8.2 


± 


1.7 








3.1 


± 


0.7 


221 


± 


16 


13 


SMMJ163658.19+410523.8 


+1.1, 


-1.2 


2.454 


45.2 ± 


5.3 


10.7 


± 


2.0 








3.4 


± 


1.1 


92 


± 


16 


14 


SMM J163704.34+410530.3 


-1.8, 


-4.4 


0.840 


21.0 ± 


4.7 


11.2 


± 


1.6 








(0.8 


± 


1.1) 


45 


± 


16 


15 


SMMJ163706.51+405313.8 


-2.3, 


+2.4 


2.374 


36.1 ± 


7.7 


11.2 


± 


2.9 








4.2 


± 


1.1 


74 


± 


23 


(16) 


SMM J163639.01+405635.9 






1.495 






5.1 


± 


1.4 








3.4 


± 


0.7 


159 


± 


27 


(17)d 


SMM J105201. 25+572445.7 






2.148 


24.1 ± 


5.5 


9.9 


± 


2.2 


4.4 


± 


1.3 


3.4 


± 


0.6 


72.1 


± 


10.2 


(18) 


SMMJ105158.02+571800.2 






2.239 






7.7 


± 


1.7 








2.9 


± 


0.7 


98.1 


± 


11.6 


(19)d,e 


SMM J105200.22+572420.2 






0.689 


15.5 ± 


5.5 


5.1 


± 


1.3 


4.0 


± 


1.3 


2.4 


± 


0.6 


57.4 


± 


13.2 


(20)d 


SMM J105227.58+572512.4 






2.142 


44.0 ± 


16.0 


4.5 


± 


1.3 


4.1 


± 


1.3 


2.8 


± 


0.5 


39.2 


± 


11.4^^ 


(21)f 


SMMJ105155.47+572312.7 






2.686 






5.7 


± 


1.7 








3.3 


± 


0.8 


46.3 


± 


10.2 



Note. — MAMBO 1200/im fluxes are from Greve et al. ( 2004), upper limits (bracketed fluxes) are from T. Greve (2006, private communication), Bolocam 
1100/im fluxes are from Laurent et al. (2006), and 850 /im SCUBA and 1.4 GHz fluxes are taken from Chapman et al. (2005). Biggs & Ivison (2006) provide 
alternative radio fluxes for many of the listed objects. The bracketed IDs indicate sources that have not been observed by the authors within the context of this 
paper, but archived data are used in the analysis. 

^SHARC-2 detection offsets are with respect to the published radio positions. 

^Both of these sources are likely contributors to the observed 1100/im flux. The Bolocam data are ignored in the analysis. 
^The alternative radio fluxes from Biggs & Ivison ( 2006) are significantly different from or more accurate than the values listed here. 
'I.SHARC-2 350 /im fluxes from Laurent et al. (2006). 

^Chapman et al. (2005): "These SMGs have double radio identifications, one lying at the tabulated redshift and a second lying at z < O.5.' 
source is assumed to be the dominant contributor. 

"These SMGs have double radio identifications, both confirmed to lie at the same redshift." 



The higher redshift 



sources, and the appropriate corrections were applied to the 
images. The final calibration is expected with high confidence 
to be more accurate than 15%, with systematic effects antici- 
pated to be less. 

The maps produced by CRUSH reveal Gaussian noise pro- 
files with an anticipated tail at positive fluxes. Hence, in 
searching for 350 /im counterparts around the published ra- 
dio positions, we calculate detection threshold s based on the 
effective number of beams A^beam ^ 1 +^Mbeam (lKovacsl2Q06h 
inside a detection area A in the smoothed image. A detection 
confidence C ^ 1 is reached at the significance level S/N, at 
which 

AT- />CX) 



1-C? 



POO 

/ ' 

J S/N 



dx. 



(1) 



As for the appropriate search radius, the probability that 
the peak of the smoothed image falls at some radius away 
from the true position is essentially the likelihood that the 
smoothed noise level at that radius is sufficient to make up for 
the deficit in the underlying signal. Therefore, for the case of 
Gaussian noise profiles, we expect the detection to lie inside 
of a 2(7 noise peak, and therefore within a maximal radius 
given by the condition S/N = 2/[l -exp(-r^^^/2cr^g^jj^)] for 
a Gaussian beam profile with abeam spread at an underlying 
signal-to-noise ratio of S/N. To this we must add the appro- 
priate pointing tolerance (at 2 a level) in quadrature, arriving 
at an expression for the detection radius around the actual po- 
sition as a function of detection significance, that is. 



r^ =4a ■ ■ 

' max ^ ^ pointing 



-^CTbeamln 



S/N 



(2) 



Note that this expression simplifies to rmax -^ 2crbeam(S/N)~^/^ 
for the case of S/N ::^ 1 and negligible cTpomting- 

The combination of S/N = 2.30 and r^ax = l^"A simulta- 
neously satisfy both constraints (eqs. Q and Q) at C = 95% 
confidence level for a abeam = 9" effective beam and apointing = 
3^'5 pointing rms. Potential candidates thus identified are sub- 
sequently verified to lie within the expected distance from 
their respective radio positions. The resulting identifications 
are summarized in Table [l] When counterparts were not 
found, the peak measurement values inside the search area 
are reported. 

The sources, collectively, are much smaller than the 
SHARC-2 beam (i.e., d < 9"), as neither (1) fitting larger 
(12^0 beams or (2) filtering extended structures produces sys- 
tematically different fluxes for the sample as a whole. Source 
extents typically < 30kpc are therefore implied. While the 
partial resolution of a few objects cannot be excluded, the 
peak fluxes of the optimally filtered images are expected to 
be generally accurate measures of the total integrated flux for 
the objects concerned. 

3. SPECTRAL ENERGY DISTRIBUTIONS OF SMGs 

We fitted the SHARC-2 350 /im a nd SCUB A 850 um fluxes, 
combined with Bol ocam 1100/im JLaurent et aL.2005.) and 
MAMBO 1200/im JGreve et alJl2004h data when available, 
with single-temperature, optically thin greybody models of 
the form Siy^ T) ex hz(i')B(u, T), where hz(u) ex v^ is an ap- 
proximation for the full emissivity term (l-expKv/vo)^]) for 
^obs "C ^0- Alternative SED models incorporating the full op- 
tical depth, or a distribution of temperatures and power-law 
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Fig. 2. — Characteristic dust temperature 7j and effective emissivity index 
f3 are highly correlated parameters of the SED fit for the typical SMG, even 
when three photometric data points allow simultaneous fitting of both Tj and 
(3 (source 13 from Table |2l shown with solid 1,2 and 3 a contours). How- 
ever, by assuming f3 = 1.5, we can accurately constrain luminosities (dotted 
contours) and temperatures, which are listed in Tables Inland |3] The derived 
properties may be scaled to other [3 values of the reader's preference using 
the relations of eq. (thick solid line). 



Wien tails, did not provide a better description of the data. 
Specifically, the flattening of the Rayleigh- Jeans slope due to 
optical depths approaching unity is not detectable with the 
:^ 10% uncertain relative calibration of the bands, while the 
Wien side of the spectra is not sampled by the observations. 
More co mplex SED models, e.g. the two-temperature model 
used bv lDunne&"E ales (2001), were not considered, since 
these require a greater number of parameters than can be de- 
termined from the few, often just two, photometric data points 
available for the typical SMG. 

SED models, whether incorporating an emissivity slope, 
full optical depth or multiple temperature components, are 
simple parametrizations of complex underlying spectral dis- 
tributions, produced by a range of dust properties inside the 
targeted galaxies. Therefore, T and (3 of the model describe 
not the physical temperature and emissivity of every single 
dust grain, but provide an effective characterization of dust 
emission in a galaxy as a whole. In this sense, the characteris- 
tic /?- values, expected in the range 1-2, reflect both the under- 
lying grain emissivities and a distribution of physical temper- 
atures within the observed objects. The derived T- values pro- 
vide an effective comparison of the characteristic dust temper- 
atures among objects with equivalent SED parametrization. 

Simultaneous fitting of T and (3 requires at least three pho- 
tometric data points (while many SMGs have only two), and 
even when permitted these are highly correlated parameters of 
the fit for the typical SMG (Fig.|2|l. Therefore, we initially as- 
sume f3= 1.5 for the SED fit, as it provides good characteriza- 
tion of actively star- forming envi ronments, b oth in accurately 
modeled Galactic clouds (Duoac et al."2003') and galaxies of 
the local universe (Dunne & Eales 2001), and is broadly con- 
sistent with laboratory measurements on carbite and silicate 
fT996) . 

via 



TABLE 2 
Far-Infrared Properties of SMGs from SHARC-2 Data 



ID 



(K) 



log Lfir 

(Lq) 



qi 



1 


43.3 ± 18.7 


12.83 ± 0.64 


2.39 ± 0.64 


1.36 


2.94 


2 


66.8 ± 26.1 


13.49 ± 0.57 


2.31 ± 0.57 


1.61 


3.69 


3 


21.3 ± 3.8 


11.95 ± 0.24 


1.91 ± 0.26 


0.50 


0.21 


4 


40.1 ± 5.8 


12.98 ± 0.19 


2.33 ± 0.21 


0.84 


1.06 


5 


39.1 ± A2 


13.03 ± 0.14 


2.39 ± 0.21 


0.69 


0.66 


6 


25.7 ± 4.7 


12.28 ± 0.18 


1.71 ± 0.18 


0.52 


0.35 


7 


31.0 ± 5.3 


12.73 ± 0.15 


2.19 ± 0.15 


0.54 


0.39 


8 


40.9 ± 8.2 


12.85 ± 0.25 


2.20 ± 0.26 


0.57 


0.46 


9 


21.1 ± 6.3 


11.95 ± 0.29 


1.16 ± 0.29 


0.31 


0.04 


10 


24.3 ± 6.8 


11.88 ± 0.35 


0.79 ± 0.35 


0.59 


0.50 


11 


53.4 ± 16.2 


13.18 ± 0.42 


2.61 ± 0.43 


1.38 


2.88 


12 


34.5 ± 4.5 


12.72 ± 0.17 


1.75 ± 0.17 


0.73 


0.77 


13 


37.3 ± 5.0 


12.93 ± 0.16 


2.31 ± 0.18 


0.77 


0.91 


14 


16.8 ± 2.3 


11.38 ± 0.15 


2.13 ± 0.22 


0.65 


0.53 


15 


31.5 ± A2 


12.71 ± 0.16 


2.23 ± 0.21 


0.65 


0.57 


(17) 


21 A ± 3.0 


12.41 ± 0.15 


2.02 ± 0.16 


0.61 


0.44 


(19) 


14.2 ± 2.2 


10.97 ± 0.21 


1.82 ± 0.23 


0.58 


0.41 


(20) 


40.6 ± 12.3 


12.87 ± 0.48 


2.76 ± 0.50 


0.99 


1.61 



Note. — A summary of the derived far-infrared properties of the observed 
SMGs. All quantities were derived using an optically thin approximation with 
13=1.5. Temperatures and luminosities may be scaled to other ^^-values using 
the relations of eq. J3j and the corresponding indices e and A liste d here. The 
qL values derived from the radio data of Biggs & Ivisoil 120061) tend to be 
higher than these by 0.06 on average. 



TABLE 3 
Properties of SMGs Incorporating Radio Data 





Td 


lOgLpiR 


log Md 






ID 


(K) 


(Lq) 


(Mq) 


e 


A 


1 


37.2 ± 3.8 


12.60 ± 0.12 


8.64 ± 0.17 


0.82 


1.09 


2 


60.3 ± 6.1 


13.34 ± 0.12 


8.22 ± 0.16 


0.86 


1.15 


3 


23.5 ± 2.5 


12.10 ± 0.15 


9.24 ± 0.17 


0.72 


0.82 


4 


37.1 ± 3.4 


12.86 ± 0.11 


8.91 ± 0.14 


0.77 


0.89 


5 


36.8 ± 3.3 


12.94 ± 0.12 


9.01 ± 0.12 


0.67 


0.65 


6 


38.6 ± 7.0 


12.66 ± 0.13 


8.61 ± 0.33 


0.94 


1.29 


7 


30.3 ± 3.9 


12.71 ± 0.09 


9.24 ± 0.25 


0.68 


0.78 


8 


39.8 ± 4.5 


12.81 ± 0.12 


8.69 ± 0.19 


0.71 


0.93 


9a 


21.1 ± 6.3 


11.95 ± 0.29 


9.36 ± 0.48 


0.31 


0.04 


10^ 


24.3 ± 6.8 


11.88 ± 0.35 


8.93 ± 0.37 


0.59 


0.50 


11 


41.5 ± 4.4 


12.82 ± 0.13 


8.60 ± 0.16 


0.80 


0.99 


12 


43.2 ± 4.7 


13.01 ± 0.10 


9.69 ± 0.17 


0.96 


1.40 


13 


34.8 ± 3.2 


12.84 ± 0.10 


9.04 ± 0.14 


0.73 


0.83 


14 


16.8 ± 1.9 


11.38 ± 0.14 


9.32 ± 0.19 


0.72 


0.74 


15 


30.7 ± 3.2 


12.68 ± 0.12 


9.18 ± 0.16 


0.70 


0.75 


(17) 


28.6 ± 2.3 


12.47 ± 0.12 


9.14 ± 0.13 


0.75 


0.85 


(19) 


16.4 ± 1.7 


11.17 ± 0.13 


9.17 ± 0.16 


0.81 


1.02 


(20) 


28.6 ± 3.5 


12.31 ± 0.17 


8.98 ± 0.15 


0.76 


1.02 



grains (lAgladze et alJ 1199(1) . The properties, thus derived 
may be scaled to the /^-values of the reader's choosing 



Note. — Quantities were derived similarly to Table|2| except that the ra- 
dio data are also incorporated with qi = 2. 14. The estimate of the dust masses 
additionally assumes ^cj(850/im) = 0.15m^kg~^ for the dust absorbtion effi- 
ciency. 

^These SMGs are identified as being radio-loud. Therefore, the radio data 
are not used in constraining the far-infrared SEDs. 

the deduced spectral indices e and A listed in the Tables |2| and 
Elas 

T, a (3-' , Lfir ^ /?" ^ , M, a /3-(4+/3)-A ^3^ 

We included 15% calibration uncertainty in addition to the 
published statistical uncertainties for all submillimeter data. 
Nevertheless the luminosities are constrained accurately be- 
cause the SHARC-2 350 /im measurement falls near the emis- 
sion peak for all of these sources (Fig.ISlK 

Flux boosting or Eddington bias ( Coppin et al 120051 ). in the 
case of less significant detections, where the selection is in 
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Fig. 3. — Optically thin greybody SED models were fitted {solid lines) to 
the various available observed fluxes and are transformed to the rest frames of 
the galaxies for comparison. Fits using ( Laurent et al. 2006 ) data also shown 
{dashed lines). All the spectra peak in the neighbourhood of the redshifted 
350 /im point, providing testimony of the pivotal importance of the SHARC- 
2 data point in constraining these models. The SEDs shown represent our 
best fit to the data for sources detected by SHARC-2, use an emissivity index 
(3 = 1.5, and incorporate radio fluxes via the far-infrared to radio correlation 
with ^L = 2.14 with an assumed 0.12 dex dispersion. 



the submillimeter, induce a small bias in the derived SEDs. 
As the bias will not apply to follow-up measurements, and 
because the exact quantification of the Eddington bias is non- 
trivial, requiring a-priori knowledge of the underlying source 
distributions, we left fluxes uncorrected for its effects. 

For the first time we are able to use multiband photometry 
to accurately determine the characteristic dust temperatures 
and far-infrared luminosities for the SMG population. Lu- 
minosities (Tables 121 and 13 are calculated analytically from 
the fitted grey body model Siy^ T), following Pe Breuck et aL 
(|200J: 

(4) 
where luminosity distances (Dl) were obtained^ for a ACDM 
cosmology with Hq = 65kms~^ Mpc~\ Qm = 0-3 and ^^a = 
0.7. The above expression provides the correct SED integral 
as long as the transition from the optically thin greybody ap- 
proximation to optically thick blackbody is above the emis- 
sion peak where the contribution to the total luminosity is 
negligible. If p ower- law Wien tails, with spectral slopes of 
-a (JBlain et al.l l2QQ3h . are assumed instead, the luminosities 
can be scaled by a constant rjpia), which is around 1.5 for 
an Arp220 type template with a = 2. More generally, in the 
range of a ~ 1 . 1-4.0 and P ~ 0-3, the values of r] are well ap- 
proximated (with an rms of 5% or 0.02 dex) by the empirical 
formula 

r]f3(a) ^ (1.44 + 0.07 /3)(a- l.OQr^^l 

^ See http://www.astro.ucla.edu/ wright/CosmoCalc.html 



Similar corrections may be derived for the case of increasing 
optical depth, with r^ as a function of hvo/kT. 

Illuminated dust masses (TablelHi were also estimated from 
the SED model S(u, T), using (De Breuck et al. 20JD, 



M,= 



S{vJ)Dl 



(1 +z)A^j(i^rest)^(^rest, Td) 

Here, the normalization for the absorption efficiency was 
assumed to be /^sso/xm = 0.15m^kg~\ representing the ex- 
trapolated average 125 /im value of 2.64 ± 0.29 m^kg~^ 
jDunne et a l."2003') from various models by assuming (3 of 
1.5. In comparison to the gas mass es derive d from^ CO mea- 
surements (Greve et al. 2005; Tacc oni et al.k2 006). we find an 
average gas-to-dust ratio of 54;^|^ (/i:850/im/0.15m^kg~^), re- 
sembling t he ratios s een in nuclear regions of local galax- 
ies by Sea auist et al.l (i2004i) . who assume a i^d^co product 
comparable to ours. The individual measurements indicate 
an intrinsic spread of c^ 40% around the mean value. The 
low ratios may be interpreted as an indication for the relative 
prevalence of dust in SMGs over the local population, which 
typically exhibit Milky Way like gas-to-dust ratios around 
120 (Stevens et al. 2005) or, alternatively, that absorption in 
SMGs is more efficient with tv^so^im ^ 0.33 m^ kg~^. 

Our estimates of the dust temperatures, and hence luminosi- 
ties, are systematically lower, than those antic ipated based 
on the 850 /im and 1.4 GHz fluxes alone by C hapman et alJ 
(2005), who overestimate these quantities by 13% and a fac- 
tor of :^ 2, respectively, when assuming the local far-infrared 
to radio correlation. This discrepancy can be fully reconciled 
if a different constant of correlation is assumed for SMGs (see 
Section|4|l. 

We confirm that the SMG population is dominated by 
extremely luminous (several times 10^^ Lq) systems with 
:^ 10^ M0 of heated dust, and characteristic 35 K dust 
temperatures typical to actively star forming ULIRGs. As 
anticipated, these objects resemble the nearby archetypal 
ULIRG, Arp 220 (with 7^ ^ 37 K, similarly obtained), except 
that they are about 7 times more luminous on average. 

3.1. Cold, Quiescent SMGs? 

In addition, there appear to be several cooler (T^ < 25 K), 
less luminous (lO^^-lO^^L©) objects, albeit with compara- 
ble dust masses (few times lO^M©) present in the population 
(sources 14 and 19, and possibly 3, 9, and 10 in Tables |2| and 
O. While these resemble the Milky Way in temperatures, and 
hence probably in star formation densities, they are tens of 
times more luminous than the Galaxy. 

This combination of extreme dust masses yet surprisingly 
low relative star formation rates, indicated by the lesser dust 
heating, allude to possibly incorrect low-redshift identifica- 
tions. Should these galaxies lie at the higher, more typical 
redshifts of the SMG population, their temperatures and lu- 
minosities would closely resemble those of the hotter pop- 
ulation. However, while lensing of distant SMGs by mas- 
sive f oreground g alaxies at the measured redshifts, is conceiv- 
able ("Blain et aUri 999^ IChanman et al J l2002). these should 
be rare, and unlikely to account for all cold SMGs. 

Alternatively, these seemingly normal type galaxies could, 
potentially, be remnants to the more remote, hotter pop- 

^ The CO fluxes from'Tacconi et al.' ^200^ were converted into gas masses 
usin g the conversion factor Xqq of O.SMo (Kkms~^ pc^ )~^ in Solomon et all 
IT997) appropriate for U LIRGs hPownes & Solomorll998.) . The results are 
therefore consistent with Greve et all I 
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Fig. 4. — Median dust temperatures vs. redshift. Data are binned for this 
paper (squares), for the Chapman et al. ( 2005) sample (circles) with tempera- 
tures estimates incorporating the radio data, and for dusty quasars (triangles). 
One-sigma error bars were estimated from the scatter inside the bins and the 
uncertainties of dust temperatures where available (thick bars). Dispersion 
within bins are also shown (thin bars), while dotted error bars indicate mea- 
surement error for single datum inside a bin. A Tj of 35 ± 3 K may be ap- 
plicable for z ~ 1.5-3.5, which, if true, can provide photometric redshifts in 
that range. However, observing frame temperatures of around 10 K appear to 
fit the SMG data somewhat better than constant 7j in the rest frame. This 
similarity of the SEDs in the observing frame may owe to selection effects. 

ulation of SMGs, once the short Hved starbursting activ- 
ity subsides. Dust could persist and remain detectable 
in the more proximate universe, provided that radiation 
from the rapidly growing black holes (Borysetal. 2005; 
[Alexander et al. 2005a) does not disperse dust entirely, e.g. 
dusty QSOs of .Benfordet al., (JJ99) and B eelen etal., (,2006 ). 
or if the dust can later reassemble to form large disk galaxies, 
like source 14, which is consistent with the large spiral disk 
galaxy seen in Hubble Sp ace Telescope images ( Almaini et al. 
|2006; Borys et al. 2006), or the cold dust mergers suggested 
by the chain of 16 /im emission seen for source 19. This con- 
nection between the hot and cold SMGs is further hinted at by 
the apparent trend of dust temperatures increasing with red- 
shift (Fig. |4ll even without the low-redshift data. We, there- 
fore, find the existence of massive cold SMGs plausible. 

4. THE FAR-IR/RADIO CORRELATION AT HIGH REDSHIFTS 

One of the most valuable uses of the new measurements 
is to test the local far-infrared to radio correlation at the red- 
shifts of SMGs (z ~ 1-3). The simple extension of the correla- 
tion into the distant universe was suggested by earlier studies 
based on a combination of 1.4 GHz radio with Infrared Space 
Observatory (ISO) 15/im ('Garrett"2002') and Spitzer 24 and 
70 /im observations ( Appleton et al. 2004), both of which take 
advantage of SED template fitting to extrapolate fluxes into 
the far-infrared regime. More recently, B eelen et aL ^2006) 
hints at the validity of the correlation in the case of a hand- 
ful of distant (z ~ 2-6) quasars with SEDs well constrained 
directly from far-infrared measurements. 

A quantitative treatment of the correlation (iHelou et alJ 
Il985h has been formulated in terms of the Infrared Astronom- 
ical Satellite (IRAS) 60 and 100 /im fluxes, expressed in the 



rest frame as 
^ = log I 



FIR 



-log 



^1 



4 GHz 



(5) 



^3.75 X 1012 Wm-V ^^"^ VWm-2Hz"V 
where the far-infrared parameter FIR = 1.26 x 
10"i'^(2.585'60/im + 5'ioo/im)Wm"2 and the fluxes S^o^m 
and ^loo/im are in Jy. The quantity essentially estimates the 
flux in a wide band centered at 80 /im and is a good tracer of 
the far-infrared luminosity for temperatures in the range of 
20-80 K and emissivity indices p -- 0-2 (Helou et al. 1981). 
Perhaps a better measure of the correlation is obtained by 
directly comparing luminosities. Therefore, we propose to 
use 



qi = log 



^FIR 



(6) 



J4.52THz]Li.4GHz 
Here the normalization frequency 4.52 THz has been derived 
for the adopted greybody model^ such that qi ^ q for T = 
40 K and (3 = 1.5, thus approximating the original definition 
of q while extending its usefulness beyond the original re- 
strictions in temperature and emissivity. Radio luminosities, 
in turn, are calculated as 

LiAGRz = 4nDlSiAGnz(l +zr-\ (7) 

which includes a bandwidth compression by (1 +z)~^ and a K- 
correction (1 +z)" to rest frame 1.4 GHz. We assume a syn- 
chrotron power law of 5" ex u'^ , with the spectral index a for 
which a typical value of 0.7 (Condon 1992) of nonthermal 
radio sources is assumed for the calculations.^ 

We confirm that the far-infrared to radio correlation ap- 
pears to hold throughout the redshift range of z ~ 1-3 (Figs.|5l 
and |6|), with the notable exceptions of sources 9 and 10, 
which likely host radio loud AGNs^^. The detections reveal 
q ^ 2.07 ± 0.09 and q^ = 2.12 ± 0.07 with intrinsic spreads 
estimated at a^ ?^ 0.21 and aq^ ?^ 0.12 around t he respective 
mean values. The alternative radio fluxes of Bi ggs & IvisonI 
uOOo) provide q values that are higher by 0.06 on average 
and scatter slightly more with aq ^ 0.30 and a^^ ^ 0.23. 

The reduced dispersion in qi vs. q and the somewhat pe- 
culiar case of the two cold lower redshift objects (14 and 
19), whose derived dust temperature of only 16 K fall out- 
side the range where FIR traces luminosity (eq. O) high- 
light the practicality of the luminosity-based ^L-measure over 
the original IRAS -hsised q. The derived mean values are 
significantly less than the locally observed median values 
of 2.75 ± 0.03 for local spirals and starbusts galaxies from 
the IRAS Faint Source Catalog JCondon & BrodericklfTool . 
2.34 ±0.01 for radi o-identifie d, flux-limited (Seoi^m > 2Jy) 
IRAS sources ( Yun et al. 2001), or :^ 2.3 for a range of nor- 
mal galaxies ( Condo n 1992) that include spirals and irreg- 
ulars, E/SO types, IRAS and radio- selected samples. The 
corresponding local scatters are also larger, with aq = 0.14, 
cr^ = 0.33, and a^ < 0.2, respectively. 

Moreover, we find the relationship to be linear, within 
measurement uncertainty, over nearly three decades of lu- 
minosities of the observed SMGs, with a best-fit FIR to ra- 
dio luminosity index, JlogLpiR/JlogLi^cHz, of 1.02 ±0.12. 

^ If SEDs include power law Wien tails instead, the normalizing frequency 
should be appropriately scaled by rip(a). Similar adjustements to the nor- 
malization can be derived for any other SED model, hence the shape of the 
particular SED model used has no effect on the correlation parameter qi oth- 
erwise 

^ The assumption of the radio spectral index a is not critical. Indices dif- 
ferent by 5a may produce a small bias in q on the order oi5q^Q.55a'vi\ the 
redshift range of z ~ 1-3 of SMGs, with an undetectable redshift dependence 
2i^5q^Gq. ^^^^^^^^_^ 



Source 9 has been characterised as an AGN byU 
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Fig. 5. — Radio to far-infrared correlation for SMGs. Galaxie s observed 
by the authors shown with squares, data from Laurent et all 12006) with cir- 
cles. Two-sigma upper limits are indicated when appropriate {arrows). Far- 
infrared luminosities were calculated exclusively based on the available sub- 
millimeter measurements, with ^5 = 1 .5. The best-fit model {solid line) reveals 
no deviation from linearity in the relationship. The deduced 2 a instrinsic 
scatters around the model are also shown {dotted lines). 

No netheless, small levels of nonli nearity , like those discussed 
by Fitt et al. (1988^; "Cox et al." ("1988); iDevereux & EalesI 
|l989) and Condon et al. ( 1991) remain possible. 

The low ^-values, the tightness of the correlation, the ob- 
served linearity, and the typically warm dust temperatures all 
point to dust heating that is dominated by high-mass (> SMq) 
stars. In the two-component model of Yun et al.^^2001^, short- 
lived (< 10^ yr) high-mass stars fuel the "hot" thermal compo- 
nent that is strongly coupled to the nonthermal radio via Type 
II supernovae, whereas the "cold" component heated by a dif- 
ferent population of lower mass stars is expected to produce 
a different correlation with a weaker coupling. They estimate 
qi = 23 and ^2 > 4 for the components. As most galaxies 
contain both components, intermediate ^-values are typically 
observed. The low q's that characterize SMGs are, therefore, 
indicative of the predominance of high-mass star formation in 
the dust heating. Moreover, the SMG result forces the revi- 
sion of the "hot" component correlation to ^1 < 2.1, or lower 
still (by about 0.16), if steeper radio spectra l slopes of a c^ 1 
are assumed for these very luminous objects J Yun et all2QQlh . 

The presence of radio-quiet AGNs would contribute to the 
total far-infrared emission, biasing towards higher values of 
q. Radio-loud AGNs, on the other hand, would produce low 
^-values, albeit to wildly varying extent, which should mani- 
fest as an increased spread in the correlation. The tightness of 
the observed correlation, therefore, confirms that radio-loud 
AGNs are rare, while the low ^-values, otherwise observed, 
indicate that the AGN contribution to the far-infared lumi- 
nosity is small when compared to heating from star forma- 
tion activity in SMGs. This is consistent with lAlexander et all 
(12005 a b), who used ultra-deep Chandra X-ray data to find 
that the AGNs, present in at least 75% of SMGs, contribute lit- 
tle to the total bolometric luminosity of their hosts. Dust heat- 
ing, therefore, is assumed to be star-formation-dominated. 

Last but not least, the low ^-values observed may par- 
tially arise from selection bias, as our initial sample would 



miss SMGs with low radio fluxes, i.e., sources with higher q- 
values, in the secondary selection required for redshift iden- 
tification. However, faint radio sources with undetectable 
submillimeter counterparts, representing a hypothetised hot 
extension to the SMG population, may be more typical 
J Chapman et al. 2004: Blain et al. 200 4). and thus selection 
bias could be reversed, missing the low ^-values instead. 

5. DUST EMISSIVITY INDEX 

It is possible to further constrain the dust properties, e.g. 
by measuring the effective emissivity index p. In princi- 
ple, this may be attempted on an individual basis, since sev- 
eral galaxies have MAMBO cGreveet al-.2004.) and/or Bolo- 
al] l2005h detections at longer wavelengths. 

In practice. 



cam JLaurent 

providing justification for a three-parameter fit. 
however, such an approach is unlikely to yield robust quan- 
tities for individual galaxies (see Fig. |2|) owing to the prox- 
imity of the observed bands and the modest significance of 
the detections. Therefore, we aimed to fit a single emissivity 
slope (3 for the entire sample, hoping to provide a better con- 
strained, ensemble-averaged, dust emi ssivity index. More- 
over, we extend our sample to Chapm an et alJ §2005) SMGs 
with multiband submillimeter data tha t were n ot observed in 
this paper. MAMBO associations (Gre veet aTl l2004) an d the 
350 /im follow-up of Bolocam and candidates ILaurent et alJ 
2006) provide five extra sources. 

We analyse the data consisting of detections only separately 
from the data including measurement peaks inside the search 
area for nondetections. While the first set may appear more 
robust, it could be more susceptible to selection biases due to 
the rigorous detection requirement. The larger set, contain- 
ing upper limits, is thus likely to be less biased of the two 
and could provide a better characterisation of the whole SMG 
population. 

We have used a modified, nested version of the downhill 
simplex method (Press et al. 1986) to perform the fit, yielding 
(3 = 2.42 ±0.71, a value that is too uncertain to distinguish 
between the expected values in the range of 1-2. However, 
the demonstrated validity of the linear far-infrared to radio 
correlation for the SMG sample allows us to incorporate the 
radio fluxes to constrain the total integrated far-infrared lumi- 
nosities and improve the quality of the fit. Here we assumed 
0.12 dex as the intrinsic spread in q^, producing the expected 
model deviation with x^ ^ I, and consi stent with the scatte r 
observed for local spirals and starbursts dCondon et alll99lh . 

To avoid radio-loud objects biasing our analysis, we ig- 
nored sources 9 and 10 due to their a priori low ^^-values. 
At the same time, we utilized a modified x deviation defini- 
tion to allow some of the sources to become "radio-loud" in 
the fit, i.e. to significantly deviate from the correlation in the 
lower qL direction. Specifically, when the derived qi value 
for a source is more than 2 a below the average ^^ of the fit, 
we use a relaxed deviation with a' = 2, more characteristic of 

radio loud objects, while keeping x^ continuous over the en- 
tire parameter space. The fraction of normal sources falsely 
identified as "radio-loud" by this method is expected to be 
around 2%, and therefore the effects of mischaracterizations 
are expected to remain limited. 

With the local median q value of 2.34 we obtain /? = 0.81 ± 
0. 19 for the detections and /?=1.00±0.19 when all the avail- 
able information is used in the analysis. Both these results 
imply an emissivity index significantly less than the usually 
quoted values around 1.5. The interpretation, however, re- 
lies entirely on the assumption that the local radio-far-infrared 
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Fig. 6. — (a) Original definition of q from Helou et al. 1 1985), and (b) the modified, luminosity-based defini tion ql as a function of redshift. Two-sigma upper 
limits are indicated as appropriate (arrows). Data from this paper are shown with squares, while those from Laurent et al.i(2006) are shown as dots. The observed 
SMG's seem to consistently fall below the locally observed mean of 2.34 ( Yun et al. 2001) for local bright IRAS sources (solid line) whose 2 a scatters are also 
indicated (dashed lines). The values were all derived assuming (3 = 1.5. The derived distribution mean values of the non-radio-loud subsample (^ > 1) are also 
shown (shaded bands). 



correlation holds unaltered for higher redshifts. 

For a more rigorous treatment, we obtained confidence con- 
tours for the likely combinations of both (3 and qi (Fig. Q, 
hence avoiding a priori assumptions on both of these quan- 
tities. The smaller set, with detection requirement imposed, 
favours somewhat lower /3- values, in combination with "nor- 
mal" qL, with the best fit located at the respective values of 
0.95 and 2.27, whereas the more inclusive data set tends to- 
wards emissivity indices in line with expectations, but with 
qi decidedly below the mean in the local universe. The best- 
fit values of the extended set are at (3 = 1.63 and ^l = 2.10. 
Accordingly, we assumed /3 = 1.5 and used the correspond- 
ing qL = 2.14 (from Fig.Q to calculate most likely estimates 
of the far-infrared emission properties, which are listed in Ta- 
bleia 

Systematic biases in the cross-calibration, between the 
SHARC-2 and the longer wavelength data points, would in- 
troduce bias into the estimated emissivity index. Fortunately, 
with expected 10% relative calibration error, the resulting bias 
is only ^/? :^ 0.08, and our conclusions on the dust emissivity 
are nearly unaffected. 

6. PHOTOMETRIC REDSHIFTS 

Obtaining accurate redshifts for distant SMGs has relied on 
optical mea surements, guided by radio or optical associations 
iChapman et al...2003, .2005.) . This approach has been very 
successful, providing accurate redshifts for nearly a hundred 
SMGs to date, but it is unlikely to be practical for compiling 
much larger samples. Moreover, the identifications involving 
radio or optical counterparts could introduce selection biases 
into the samples. For both reasons a "holy grail" of large sub- 
millimeter surveys has been to obtain redshifts directly from 
the radio, submillimeter, and mid- to far-infrared (e.g. Spitzer) 
if possible. 

The simple st types of redshift estimators, such as 
Canlli&Yunl (ITQ9Q. 2000a b), Wiklind OQQD and 
Laurent et al.l (l2006h . assume that all SMGs have intrin- 



sically similar rest- frame temperatures. In essence, redshifts 
are calculated as {\+z)^ Td/Tohs (Blain 1999) for SEDs 
characterized by Jobs in the observing frame and an as- 
sumed representative dust temperature T^ for all objects. 
The methods only differ in how Tobs is estimated from the 
various available flux data. Only independent, and accurate. 




Fig. 7. — Likelihood contours (at 68%, 90%, and 95% confidence) for the 
dust emissivity [3 and the radio to far-infrared correlation constant qi, using 
only the available detections (dashed curves) and all information, including 
nondetections (solid curves). The best-fit locii of both fits are indicated by 
crosses. The more inclusive set is expected to be a more reliable indicator of 
SMG properties, as it is less affected by the selection biases that are intro- 
duced by requiring detections in all bands. A similar fit to Arp 220 is also 
indicated (black dot). 
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Fig. 8. — Luminosity-temperature (L-T) relation. One-sigma likelihood contours for greybody models fitted to submillimeter data (this paper with solid curves, 
Laurent et al. [2006] as dotted curves). Contours are also shown for objects not detected by SHARC-2 (dashed). The dusty quasars (triangles) of Benford et al.] 
fl999) and Beelen et al. ( 2006), and some well-known objects with accurately constrained SEDs are also shown. The plot highlights the difficulty of measuring 
the underlying temperature dependence of luminosity, as the two parameters are correlated; most data points are heavily clustered and clearly affected by selection 
effects. When such a relation is established, however, it can dispute the idea of luminosity-based photometric redshifts. The top graph shows the apparent L-T 
relation and the 850 /im selection functions (^Sgso/^m > 5mJy) at redshifts of 0.5, 1, 2, 3, and 4. (a) The observed relation is evidently biased by selection, (b) 
A bias- and evolution-free version shows the same relation after the observed redshift evolution (including selection bias) has been accounted for, effectively 
"moving" luminosities into the present era. The case of nearby Arp 220 clearly supports the case of explixit redshift evolution of the L-T relationship. 



measurements of actual dust temperatures for SMGs with 
known redshifts can test the vaHdity of this approach. We 
offer this test for the first time. 

A glance at Tables|2|and[3lreveals that dust temperatures are 
distributed around 35 K, with some 80% of the data between 
25-45 K. The same conclusion is suggested by the larger 
[chapman et al. (2005) sample, with the dust temperatures es- 
timated from the SCUBA 850 /im and radio data using the ap- 
propriate far-infrared to radio correlation as discussed above. 
This majority population of SMGs is consistent, within mea- 
surement uncertainty, with a single representative temperature 
of 34.6 ± 1 .4 K in the rest frame, with an estimated 3 K intrin- 
sic dispersion around the mean, lending some credibility to 
photometric redshifts derived as 



(1+z)^ 



(34.6±3.0)K /1.5 



0.71 



in the redshift the range of z ~ 1 .5-3.5. However, the observa- 
tions suggest that these photometric redshift indicators could 
be wholly inappropriate for as many as 20% of all SMGs, 
unless some of these outliers have incorrect redshifts. Curi- 
ously, the very low g^ -valu es measured for sources 14 and 19 
hint at "hot" type galaxies JYun et alJl2001h . in apparent con- 
tradiction with the cool dust temperatures that are implicated 
at their low redshifts. This mismatch would be reconciled if 
the temperatures were typical at the photometric z- values of 
2.96!^-^J and 1.95!^-?^ respectively. 



-0.43 ' 



While quasars (Benford et al. 1999; Be elen et aH2006h too 
can be characterised by rest- frame SED templates, with 7j ^ 
48 ± 8K^^ more or less uniformly across z ~ 2-7, they are 
typically hotter than SMGs, perhaps due to additional dust 
heating by an AGN. However, quasars, selected from larger 

006f) were adjusted for f3 = 1.5 for 



^^ Temperatures from lBeelen t 
consistency. 



volumes than SMGs, could provide the rare extension of the 
SMG population into ever more extreme luminosities and hot- 
ter dust temperatures. If true, the different temperatures that 
characterize quasars versus the bulk of SMGs invalidates the 
single representative temperature assumption for all SMGs, 
already under fire from a minority population of cold SMGs, 
and therefore may oppose rather than support the applicability 
of temperature-based photometric redshifts. 

Alternatively, the entire SMGs sample, including the low- 
redshift data, is better fit by temperature evolution of Td ex 
(1 +z), producing similar temperatures in the observed frame 
(Jobs — 10 K), with an average deviation of |x| = 1.40 versus 
2.67 produced by the representative rest- frame temperature 
model when assuming a 10% intrinsic scatter in T^. Hence, 
the temperature range that fits most SMGs could arise from 
the similar fractional dispersion of the redshift distribution, 
i.e., objects selected near the median redshift of 2.3 tend to 
have temperatures of about 35 K. Thus, temperature-based 
photometric redshifts may not provide additional information 
beyond a reflection of the median z of the sample. 

A similarity of SEDs in the observed frame, yielding a 
nearly constant Tobs and thus 7^ ~ (1+z), could arise sim- 
ply from selection effects that result from the minimum flux 
requirement for detections (Fig. [8|), and the small dynamic 
range (of a factor :^ 3) of flux densities that characterise cur- 
rent SMG samples. As such, selection effects alone could ren- 
der far-infrared and radio based photometric redshift indica- 
tors wholly unsuitable for unbiased flux- and volume-limited 
samples of SMGs. 

6.1. Luminosity -Temperature Relation 

An alternative, luminosity-based approach exploits the hy- 
pothetical relationship between dust temperatures and the 
luminosities that they fuel, to derive photometric redshifts 
by comparing the rest-frame relationship to observed tem- 
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perature indicators and fluxes. Like other scaling relations 
that characterize galaxies, the L-T relationship would, pre- 
sumably, be a genuine property of galaxies, reflecting the 
physics of galactic structure and dynamics. While the idea 
has been drafted by Yun & CariUi (2002), Blain et al. (2003), 
and Aretxagaetal. (2005), the particular implementation of 
such a method was left awaiting the determination of the ac- 
tual form of such an L-T relation. 

We find that the far-infrared luminosity is well approxi- 
mated by a power law of the form 



^FIR ■ 



■.C^T]x{\+zr 



(8) 



for some exponents 7 and ji. The expression incorporates the 
possibility of explicit redshift evolution, approximated by the 
term (1 +z)^ chosen for convenience. This independent red- 
shift term is expected to absorb selection effects and allow the 
determination of the underlying temperature dependence. 

The possible combinations of Lfir and T^ pairs are corre- 
lated (Fig.O, which adds difficulty to constraining the under- 
lying relationship. With this warning in mind, our best esti- 
mates of the exponents are 7 = 2.76 ±0.29 and /i = 1 .74 ± 0.38 
for aU data points and 7 = 2.61 ±0.30 with fi = 1.25 ±0.48 ex- 
cluding the lower redshift points. The corresponding normal- 
izations, £0, are 3.4 x 10^ L© (±0.33 dex) and 1.1 x 10^ L© 
(±0.45 dex), respectively. The fits were obtained via a two- 
variable linear regression of the logarithms with error bars 
(Tables on both the logL and log T values. As most galaxies 
lie in a narrow temperature and luminosity range, we point 
out that the derived temperature exponents are constrained 
by relatively few points, some or all of which could be pe- 
culiar, and therefore the true uncertainties of 7 and /i are 
likely underestimated. Nonetheles s, these scaling models suc- 
cessfully apply to the quasars of iBenford et al.l (119991) and 
^eelen et al. (2006), and the local ULIRGs Arp 220 and Mrk 
231. If we consolidate our estimates with the inclusion of 
these objects, we obtain 7 = 2.82 ±0.29, /i = 1.59 ±0.18, and 
Co = 3.6 X 10^ L0 (±0.46 dex). This relationship, derived for 
SMGs, quasars, and ULIRGs, could hint at a possible connec- 
tion among these types of objects; however, it may not apply 
to the less active galaxies of the local universe (Fig.[8ll. 

Assuming the exponents and scaling are properly deter- 
mined, we may attempt to derive a photometric redshift in- 
dicator by comparing the expressions for the far-infrared lu- 
minosity from equations ^ and ^ and from the definitions 
in equations ^ and (Q. In order to facilitate the derivation 
of analytic solutions, we approximate the luminosity distance 
by Di ^ Z)o(l -^z)^, the effect of which is to simplify all red- 
shift dependence to powers of (1 +z). This simplification is 
reasonable for the study of SMGs (leading to an rms Sz/z of 
:^ 6%) in the redshift range of 0.5-4 with Dq = 1.63 Gpc and 
J = 2.07 ± 0.03. We find that the far-infrared luminosity is 
then proportional to 



> 1.4 GHz 



(1+Z> 



25+1- 



^TUl+z) 



l+fJ' 



(X 



^ T 



4+/3 
obs 



^hu/kT _ I 



(1+Z) 



25 



(9) 

Substituting the quantities assumed or determined so far, 
we find that all three expressions possess similar dependences 
on (1+z), with exponents of 4.44 ±0.12, 4.41 ±0.34, and 
4. 14 ± 0.06, respectively. As the exponents are indistinguish- 
able within the measurement uncertainties involved, we con- 
clude that, in the redshift range of z ~ 0.5-4 of the approxima- 
tion, it is very unlikely that a meaningful measure of redshift 
can be deduced in this way — bad news for luminosity-based 



photometric redshift estimation. While it may simply be an 
unfortunate coincidence, it is more likely that the inability to 
determine redshifts in this manner is, once again, the product 
of the sample selection yielding candidates with similar ob- 
served SEDs. With these being nearly identical in the observ- 
ing frame, photometric redshifts of any kind remain elusive. 

We can, nevertheless, derive a true luminosity-based pho- 
tometric redshift indicator in the low-redshift limit z <C 1, 
where the luminosity distance takes the asymptotic form 
Dl -^ cz/Hq. In this limit. 



1 /- ^0 

IGpc — 

c 



^077 

^IGpc 



1/2 



where Li opc is what the luminosity would be if the galaxy 
were at a distance of 1 Gpc. While this is irrelevant to the 
SMG sample, it provides a useful expression for the study of 
the L-T relation in the nearby (z <C 1) universe, when values 
of z are already known. 

Comparing the L-T relationship (eq. |8|) with the expres- 
sion for the integrated far-infrared luminosity (eq. l4l) we 
infer that Md ex T^^ +/^~^)q +zY, i.e., at any given redshift 

Td ex m]^^^^^~^\ that is Td ex m/-^^^^-^\ which is expected to 
be free of selection bias. The dust heating is predominantly 
fueled by the luminosity of massive (> 8M0) stars, and so 
we find that the number of high-mass stars A^hm ex 7^ , assum- 
ing that the high-mass IMF and stellar content of all SMGs 
is similar. We may estimate the high-mass star formation 
efficiency in SMGs by t/hm oc Nuu/Md ex M^^^ +p-^).^ ^^^ 

7/HM ex M^ 

Surprisingly, this implies that more dusty galaxies produce, 
and thus contain, fewer massive stars relative to dust mass. 
Perhaps this is due to more dusty environments favoring lower 
mass star formation. This conclusion is consistent with the 
evolution of q from lower values in SMGs to higher values 
in nearby dusty disk galaxies ( Yu n et al. 2001) and with the 
dominance of high-mass stars in SMGs giving way to more 
radio-quiet low-mass stars heating dust in the local universe 
CBlaineta]_199 9b: Baugh et al. 2005). 

7. SCALING RELATIONS 

After cancelling the nearly identical redshift terms appear- 
ing in the luminosity scaling relations (eq. |9|), we obtain a 
useful set of expressions interrelating the observed quantities, 
^L4GHz, Su and Tobs in the far-infrared. For distant SMGs, in 
the redshift range z ~ 0.5-4, we find 



^L 



4 GHz 



ocr 



obs' 






{e"'' 



/kT _ 



The appropriate constants of proportionality are easily derived 
from the relevant expressions (eqs. ID, Q, |6| and |8l) ap- 
plied at the median redshift of the sample. Figure|9ltestifies to 
the predictive power of these relations. Excluding the radio- 
loud points, which clearly lie away from the far-infrared to 
radio correlation, the model is consistent with observations 
within the measurement uncertainties. We emphasize, that 
these relations are not fitted to the data, but rather a direct 
consequence of the radio to far-infrared correlation, and the 
deduced L-T relation. ^^ Therefore, the excellent agreement 
with data is remarkable. 

^^ A pure L-T relation without the evolution term yields a poor model of 
the observed relations which therefore directly support the inclusion of the 
explicit redshift evolution in eq. |8|. 
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Fig. 9. — Relations between observed quantities. («) Radio and far-infrared 
fluxes, (b) radio fluxes, and (c) far-infrared fluxes vs. observed temperatures. 
Results from this paper are plotted as squares with upper limits indicated as 
appropriate {arrows), Laurent et al. (2006) data are shown as circles. The 
model deriving from the radio-FIR correlation and the observed L-T relation 
is also shown {solid curve). The data show excellent agreement with predic- 
tions within the measurement uncertainties involved. A single radio or sub- 
millimeter measurement may, therefore, suffice to predict the other quantities. 
Thus, these relationship can offer a useful tool for all radio or submillimeter 
surveys. 



8. CONCLUSIONS 

New 350 /im data lead to the first direct measures of the 
characteristic dust temperatures and far-infrared luminosities 
of SMGs and the following valuable conclusions. 

1 . The linear radio to far-infared correlation remains valid 
out to redshifts of z ~ 1-3 for SMGs, with the excep- 
tion of rare radio-loud objects. The power-law index in 
the correlation is 1 .02 ± 0. 12, and tighter than observed 
locally, with an intrinsic dispersion of only 0.12dex. 

2. Either the far-infrared-radio correlation constant q is 
lower than locally measured {q^ ~ 2.14 ± 0.07) or the 
effective dust emissivity index is less, with /? ^ 1 con- 
sistent with observations. The lower (7- value indicates 
that dust heating is dominated by high-mass stars and 
that any AGN contribution is minor in comparison. 

3. Compared with low-redshift galaxies, SMGs are char- 
acterized either by low gas-to-dust ratios, around 
54;^|i(Ai:850/im/0-15m^kg~^), indicating dust-rich envi- 
ronments, or by efficient dust absorption of A^gsop-^ -^ 
0.33m2kg-i. 



J/im 



4. Far-infrared- and radio-based photometric redshifts 
might be appropriate for up to 80% of SMGs in the 
redshift range of z ~ 1 .5-3.5, with a rest- frame temper- 
ature assumption of 34.6 ± 3.0K(1.5//3)^-^^ However, 
photometric redshift indicators may not be appropriate 
for an unbiased SMG sample, as existing selection ef- 
fects tend to favor measuring SEDs that are similar in 
the observed frame. 

5. We deduce an L-T relationship that we argue to be a 
property of star- forming galaxies. When recast as mass 
scaling, we find, that ISM heating scales inversely with 
the quantity of dust in the galaxies i.e., the galaxies with 
more dust are less active, possibly due to the formation 
of more low-mass stars being favoured in dustier envi- 
ronments. 

6. From the observed L-T relation, possibly biased by se- 
lection, and the radio to far-infrared correlation, we 
derive scaling relations among the observed quantities 
^1.4 GHz, ^v in the submillimeter or far-infrared, and the 
observing frame dust temperature Jobs, applicable to the 
redshift range of z ~ 0.5-4. A determination of one 
of these quantities may sufficiently characterise the ob- 
served far-infrared properties of SMGs. 
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